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A driven resonator with focusing Kerr nonlinearity shows stable localized structures in a region far 
from modulational instability. The stabilization mechanism is based on front interaction in bistable 
regime with an inhomogeneous injected field. The experimental setup consist of a focusing Kerr 
resonator filled with a liquid crystal and operates in negative optical diffraction regime. Engineering 
diffraction is an appealing challenging topic in relation with left-handed materials. We solve the 
visible range of current left-handed materials to show that localized structures in a focusing Kerr 
Fabry-Perot cavity submitted to negative optical feedback are propagating fronts between two stable 
states. We evidenced analytically, numerically, and experimentally that these fronts stop due to the 
spatial inhomogeneity induced by the laser Gaussian forcing, which changes spatially the relativity 
stability between the connected states. 

PACS numbers: 


I. INTRODUCTION 

Localized structures (LS’s) often called cavity solitons in dissipative media have been observed in various fields of 
nonlinear science (see last overview on this issue [1-3]). Localized structures consist of isolated or randomly distributed 
peaks surrounded by regions in the homogeneous steady state. Currently they attract growing interest in optics due 
to potential applications for all-optical control of light, optical storage, and information processing. They appear in 
an optical resonator containing a third order nonlinear media such as liquid crystals, and driven coherently by an 
injected beam in a positive diffraction regime and close to modulational instability [4, 5]. The existence of localized 
structures and localized patterns, due to the occurrence of a modulational instability, has been abundantly discussed 
and is by now fairly well understood [6]. In this case, LS’s appears in the subcritical modulational instability regime 
where there is a coexistence between the homogeneous steady state and the spatially periodic pattern. Localized 
structures consist of isolated or randomly distributed spots surrounded by regions in the uniform state. They may 
consist of peaks or dips embedded in the homogeneous background. 

However, localized structures could be formed in modulationally stable regime [7]. In this case, heterogeneous 
initial conditions usually caused by the inherent fluctuations generate spatial domains, which are separated by their 
respective interfaces often called front solutions or interfaces or domain walls [8] . Interfaces between these metastable 
states appear in the form of propagating fronts and give rise to a rich spatiotemporal dynamics [9, 10]. Erom the 
point of view of dynamical system theory at least in one spatial dimension a front is a nonlinear solution that is 
identified in the comoving frame system as a heteroclinic orbit linking two spatially extended states [11, 12]. The 
dynamics of the interface depends on the nature of the states that are connected. In the case of a front connecting a 
stable and an unstable state, it is called as Eisher-Kolmogorov-Petrosvky-Piskunov (EKPP) front [13-16]. One of the 
characteristic features of these fronts is that the speed is not unique, nonetheless determined by the initial conditions. 
When the initial condition is bounded, after a transient, two counter propagative fronts with the minimum asymptotic 
speed emerge [13, 16]. In case that the nonlinearities are weak, this minimum speed is determined by the linear or 
marginal-stability criterion and fronts are usually referred to as pulled [16]. In the opposite case, the asymptotic 
speed can only be determined by nonlinear methods and fronts are referred to as pushed [16]. The above scenario 
changes completely for a front connecting two stable states. In the case of two uniform states, a gradient system 
tends to develop the most stable state, in order to minimize its energy, so that the front always propagates toward the 
most energetically favored state [17]. It exists only as one point in parameter space for which the front is motionless, 
which is usually called the Maxwell point, and is the point for which the two states have exactly the same energy for 
variational systems [18]. 

In this chapter, we present analytical, numerical and experimental investigation of of front interaction in bistable 
regime with an inhomogeneous injected beam. Ear from any modulational instability, the inhomogeneous injected 
beam in the form of Gaussian, can lead to the stabilization of LS’s. This completes our previous communication of 
the stability of LS’s [19]. We consider an experimental setup which consists of a focusing Kerr resonator filled with a 
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Figure 1 : Inhomogeneous Kerr cavity with negative diffraction: (a) Experimental setup, (b) Schematic representation of the 
physical cavity and the equivalent cavity, with negative cavity length. 01 optical isolator; BS beam splitter; EOM electro-optical 
modulator; GP Gian polarizer; Ci and Ci cylindrical lenses; LG liquid crystal slice; Di and D2 photodetectors; MM motorized 
mirror; Li and L2 lenses of focal length /: p and s are the polarized components of the pump (solid line) and probe (dashed 
line) beams respectively; Cnf and Cff near field and far field cameras; Mi and M2 are the real cavity mirrors but the optical 
Perot-Eabry cavity is delimited by Mi and M'2 mirrors and its length is d. 


liquid crystal and operates in negative optical diffraction regime. Engineering diffraction is an appealing challenging 
topic in relation with left-handed materials. We solve the visible range of current left-handed materials [20-22, 24] to 
show that LS’s in a focusing Kerr Fabry-Perot cavity submitted to negative optical feedback are propagating fronts 
between two stable states. Moreover, we have evidenced that these fronts stop due to the spatial inhomogeneity 
induced by the laser Gaussian forcing, which changes spatially the relativity stability between the connected states. 

This chapter is organized as follows. After an introduction we present the experimental setup and experimental 
observation in Sec. 2. The analytical and numerical investigations are presented in Sec 3. The special case with a 
null diffraction is developed in Sec. 4. We conclude in Sec. 5. 
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Figure 2: Transfer function of the optical cavity for : pump beam (blue line) and probe beam (red line), (p linear phase shift 
of the cavity. Sd cavity length variation. 


II. EXPERIMENTAL 
A. Set-up 

The experiments have been carried out using a nonlinear Kerr slice medium inserted in an optical Fabry-Perot 
resonator. The Kerr focusing medium is a 50-/im-thick layer of E7 nematic liquid crystal homeotropically anchored. 
Two plane mirrors Mi and M 2 define the physical cavity but the optical one is delimited by Mi and M '2 (which is the 
image of M 2 through the 4/ lens arrangement) and its optical length is d [Fig. 1]. The intensity reflexion coefficients 
of Ml and M 2 mirrors are Ri = 81.4 and R 2 = 81.8 respectively so that the cavity finesse is estimated to 15. The 
experimental recording of the Airy function (blue curve) on Fig. 2 gives a finesse of 11.6 indicating the presence of 
supplementary losses due to nonlinear medium and lenses transmission coefficients. 

The optical cavity length d may be tuned from positive to negative values [positive on Fig. 1(a) and negative on 
Fig. 1(b)]. Thus, for negative optical cavity path (d < 0), a beam propagating along this path experiences negative 
diffraction. Together with the positive Kerr index, the experimental setup is then equivalent to a Kerr cavity that 
would have a positive optical distance but negative Kerr index (the rja product sign in Ref. [26] that defines the type 
of transverse instabilities remains the same). However the physical mechanisms of negative refraction and negative 
diffraction are different. Thus, this intra-cavity geometrical lens arrangement allows for achieving an equivalent 
left-handed Kerr material in the visible range. It also allows to continuously tune the diffraction from positive to 
negative. 

The laser source used is a single mode frequency doubled Ad^+:YV04 laser (Aq = 532 nm) at bW. We have split 
the laser beam into two parts, a pump beam containing 95% of the initial power and a probe beam with 5% of the 
initial power. The pump beam power is controlled by an electro-optical modulator associated with a Gian polarizer. 
Then it is shaped by means of two cylindrical telescopes. The resulting beam size (^ 200 /im x 2800 jam) gives a 
"cigar” transverse shape such that only one spot can develop in one of the two directions. This beam is injected inside 
the cavity may then be considered as one-dimensional. The probe beam polarization is rotated to 90° to prevent 
interference phenomenon inside the cavity between the two beams (see Fig. (1)) The pump beam propagating in the 
forward direction is monitored at the output of M 2 mirror. We record the near field on Cnf and the far field on 
CpF- The cavity detuning is an important parameter, because it determines the solution type (mono/bistable) of the 
cavity and also the energy quantity inside the cavity, as shown on the cavity transfer function on Fig. 2. For these 
reasons, we perform an active stabilization of the cavity phase shift, based on Coen’s work and coworkers [25]. We 
stabilize the cavity length around a reference phase on the beam probe {pref on Fig. 2), then we measure the phase 
shift between the probe and the pump beams to have the real cavity detuning p, as presented on Fig. (2). Our setup 
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Figure 3: Temporal evolution of cavity transfer function correlated to cavity phase shift: (a) without stabilization; (b) with 
stabilization, d = 5 mm, / = 37 Wjerr?. 

allow to maintain the detuning long enough to realize experiments, it means to have a fix value ±7r/19 during many 
minutes. Without the stabilization, the cavity detuning can evolve quickly and not permit us to do experiments. The 
typical detuning evolution without and with an active stabilization are presented on Fig. (3). 


B. Experimental observations 

As the input power is suddenly increased to the upper bistability response branch (t = 0 s), the central part of 
transmitted intensity profile suddenly jumps after some latency time [t 37 s on Fig. 4(c)] to a higher value and 
invades the surroundings towards the external regions where the field is less intense before stopping its propagation. 
Finally, the fronts locks to give a localized light state [Fig. 4(b)]. Changing the waist w of the Gaussian forcing or 
its intensity within the bistability region allows to tune the distance between the bounded fronts and so the localized 
state extension. 

A transverse cut of the transmitted intensity profile is depicted on Fig. 4(a) obtained in the initial and the final 
observation period. This figure emphasizes the coexistence of different states in the same region of parameters. In 
addition, a state emerges from the other because of the inherent fluctuations of the system. Figure 4(c) shows this 
phenomenon starting at t « 37 s. From this instant, the system exhibits two counter-propagating fronts, which will 
become asymptotically motionless [see Fig. 4(c)]. The motionless front is observed at the location XQ{expt) = ±0.17u;. 
At this location, the input intensity is only 3 % lower than at the center of the input beam. Thus, we have observed 
experimentally coexistence between two inhomogeneous states, the noise induces a pair of fronts among these states, 
which initially are counter-propagate and asymptotically they stop. 


III. ANALYTICAL RESULTS 
A. Model equation 

The dynamics of the single longitudinal mode of the bistable system which consist of a Fabry-Perot cavity filled with 
liquid crystal Kerr like medium and driven by a coherent plane-wave steady can be described by the simple partial 
differential equation (the LL model [26]) in which we incorporate an an inhomogeneous injected field. This mean field 
approach model is valid under the following approximations: the cavity possess a high Fresnel number z.e., large- 
aspect-ratio system and we assume that the cavity is much shorter than the diffraction and the nonlinearity spatial 
scales; and for the sake of simplicity, we assume a single longitudinal mode operation. Under these assumptions the 
space-time evolution of the intracavity field is described by the following partial differential equation 

— = EM - (1 + iA)E + i\E\^E - * |a| — + (x, t) (1) 

which includes the effect of diffraction, which is proportional to a; £" is the normalized slowly-varying envelope of 
the electric field, A is the detuning parameter, Ein is the input field assumed to be real, positive and spatially 
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Figure 4: Experimental front propagation in negative diffractive inhomogeneous Kerr cavity, (a) Transverse cross-section of 
the initial (final) average localized structure in dashed black line (in continue gray line), (c) Spatiotemporal response to a step 
function of the input intensity from the lower to the upper branch of the bistable cycle, (b) cross section of the localized state. 
Light region account for high intensity of the light. Iq — 433 d = —5mm, (f = —0.6 rad, Wx = 1400 /xm, Wy = 100 /xm, 

= 81.8%; R2 = 81.4%. 


inhomogeneous. The negative diffraction coefficient is \a\. Note that the above model has been derived for a cavity 
filled with left-handed material operating in negative diffraction regime [21]. 5 scales the noise amplitude and ^(x,t) 
are Gaussian stochastic processes of zero mean and delta correlation introduced to model thermal noise [27]. 


B. Derivation of F-KPP equation 

The following development is realized for the deterministic case of Eq. 1 (5 = 0). The homogeneous steady states 
of Eq. (1) are solutions of Ei = [1 i{A — \Es\‘^)]Es. The response curve involving the intracavity intensity \Es\‘^ 
as a function of the input intensity \Ei\‘^ is monostable for A < ^/3 and exhibits a bistable behavior when the 
detuning A > ^/3. For A = ^/3, we obtain the critical point associated with bistability where the output versus input 
characteristics have an infinite slope. At the critical point, the coordinate of the intracavity are Ec = Uc + ivc with 
Uc = I\f2 and Vc — —1/3^/^\/2, and the injected field amplitude is = 2\/2/3^/^. 

The analytical investigation of fronts dynamics connecting two-homogeneous steady states in the framework of the 
Lugiato-Lefever model Eq. (1) is far from the scope of the present chapter. In this section we perform a derivation 
of a simple bistable model with inhomogeneous injection to study analytically the dynamics of the front connecting 
the two-homogeneous steady states. To do that, we introduce a small parameter that measures the distance from the 
critical point C ^ we express the cavity detuning in the form 

A = Ac(l + (2) 

where a is a quantity of order one. Then, we decompose the envelope of the electric field into its real and imaginary 
parts: E = xi-\- ix 2 and we introduce a new space and time scales as (x, t) = 3i(x/ (v^)^ • The injected field 

can be expanded as 


E — E^ 
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(3) 
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Let {u^v) = {xi^X 2 ) — {uc^Vc) be the deviations of the real, the imaginary parts of intracavity field with respect to 
the values of these quantities at the critical point. Inserting these expansion into the LL model and using the above 
scalings, we obtain 


du 

dt 

dv 

m 


+ — iv ^ + — 


Aa 


v?v 


1 d^v 
v/3 dx‘^ ’ 


-2{u + + C 



^2 

+ —{6Iu — ?>v? 


uv‘ 


■) + V3 


dx‘^ 


( 4 ) 

( 5 ) 


Our aim is to seek solutions of Eqs. (4,5) in the neighborhood of the critical point associated with the optical 
bistability. To this end, we expand the cavity field and the injected field as 


{u, v) = C[(wo, t’o) + C(wi, + C^(W 2 , V2) + C^(W3, V3) + ...] 


(6) 


Inserting these expansions and taking into account of Eqs. (2,3) into the LL model and using the above scalings, we 
then obtain a hierarchy of linear problems for the unknown functions. At the first order in C, we find uq = —vq. At 
the second, we have vi = —ui + 3//2 — Uq. Einally, at the third order, we get a simple bistable model 


du 

dt 


f{u) 


d^u 

dx‘^ 


( 7 ) 


where u{x^t) = ^3/ {2a)uo is a scalar field that accounts for the real part of the envelope E, 

f{u) = T] U — 


( 8 ) 


and T] = 4^2C^/3 controls the relative stability between the equilibria. Note that ^2 is proportional to the pumping 
Ein. Hence, if the pumping is inhomogeneous then the parameter 77 is also inhomogeneous. Eor a Gaussian pumping, 
we consider 


rj{x) = 7 ) + r]oe , (9) 

where 770 accounts for the strength of the spatial pumping beam and w is the width of the Gaussian. Eor r]{x) = 0 
both states are symmetric corresponding to the Maxwell point, where a front between these states is motionless. 

To perform analytical developments, we use further approximation by taking into account only the first order 
development close to the center of the optical pumping where the stress is maximum, i.e., 

ri{x) « + 7?o (1 - (x/wf) (10) 

close to the Maxwell point, one can consider the following ansatz for the front solution u{x,t) = 
tanh [(x — xo(t))/\/2] + iL, where xq is the front position and H accounts for small corrections. We substitute 
the above ansatz for u{x^t) in Eq. (7), linearizing in H and imposing the solvability condition, we obtain the kinetic 
equation for the evolution of the front position 


Xq = 


-3^2 

2 



( 11 ) 


this equation takes into account of corrections imputable to inhomogeneities of the injected signal. The stationary 
solutions of this equation reads 


Xq 




the term 7r^/(67n) is negligible for w large. In this case, we have 


( 12 ) 



Xq = Ew 


(13) 
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Figure 5: Numerical study of front pinning with a parabolic forcing, (a) Spatiotemporal diagram; (b) Transverse profile of the 
front at t = 4000; (c) Relative front position xo{t) (for x > 0); (d) front velocity. Blue diamonds are the numerical points and 
the red curves are the analytical predictions .770 = 0.6, fj = 0.3, w = 350. 


The ordinary differential equation (11) admits an exact solution that leads to the following trajectory of the front 


xo(t) = ±a tanh (6(t — to)), (14) 

with a, b and to are coefficients depending on 770 5 V The equilibrium position of the front xooo can be inferred 

from expression Eq. (14) for t ^ oc as xooo = y/l — 77 / 770 . This equilibrium positions could be obtained directly 
from Eq. (7), assuming rj{xooo) = 0 which is the condition for motionless front (Maxwell point). Extending this last 
property to the initial Gaussian forcing, we get 


Xooo = 



(15) 


At leading order, Eq. (15) recover again the previous expression of xqoo for the parabolic approximation of the Gaussian 
profile. 


C. Numerical results 

We conduct numerical simulations of the imperfect pitchfork bifurcation, Eq. (7-8), with a parabolic spatial injection 
to compare with the analytical predictions. Initially, we observe the front propagation, as the case with a plane wave. 
Then, the front speed decreases to become zero, around t = 1000, as presented on [Eig. 5(a,c,d)]. At t = 4000, the 
structure is pinned by the pump parabolic profile and presents a bistable structure with a parabolic dependance on 
the two states [Eig. 5(b)]. We plot on the same graph [Eig. 5(c,d)], the analytical predictions of the core front position 
and its speed with the numerical simulations. We have an excellent agreement between both. Eor the Gaussian case, 
which is not presented in this document, the hyperbolic tangent allows always to reproduce the front trajectory. We 
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Figure 6: Front propagation in negative diffractive inhomogeneous Kerr cavity. (a,b) Transverse cross-section of the initial 
(final) average localized structure in dashed black line (in continue gray line). (c,d) Spatiotemporal response to a step function 
of the input intensity from the lower to the upper branch of the bistable cycle. (e,f) Experimental and numerical cross section 
of the localized state. Light region account for high intensity of the light.(a,c,e) Experiments /o = 433 d = —5 mm, 

(p = —0.6 rad, Wx = 1400/xm, Wy = 100/xm, = 81.8%; R 2 = 81.4%. (b,d) Numerical simulation of LL model with Eq = 1.9, 
A = 3.0, a = 0.001, Wx = 1400 /xm, Wy = 100 /xm, £ = 0.4. 


need only to adjust the parameters a and b from Eq. (14) to have a good agreement between the predictions and the 
numerical simulations. 

For the LL model, Eq. (1), we focus our numerical investigations with a bi-dimensional Gaussian spatial injection 
and by taking into account the experimental parameters are shown in Fig. 6. 

We perform numerical simulations with a asymmetric Gaussian forcing with cigar shape {wx ^ Wy) 


^in y) 


Eoe 




(16) 


Furthermore, we take account the inherent fluctuations of the system by the stochastic Lugiato-Lefever model described 
by Eq. 1, 5 0. We use a stochastic Runge-Kutta solver of the order of 2 with additive noise [28]. In this latter case, 

the temporal step (At) is equal to 0.01. Numerical simulations with these ingredients show a quite good agreement 
with the experimental observations (see Fig. 4 and Fig. 6). Hence, the analytical expression (14) can be used to 
figure out and to characterize the experimental front dynamics. Figure (7) depicts the experimental and numerical 
temporal evolutions of the front position. It clearly evidences that the expression of Eq. (14) reproduces well the front 
dynamics. Therefore, the effect of a spatial forcing on front propagation is to induce the front moves and stops on an 
asymptotic position, satisfying a hyperbolic tangent trajectory. 


IV. ZERO DIFFRACTION CAVITY 

We have explored an interesting limit case, where the diffraction is zero. This limit case delimits the positive 
diffraction case where exists soliton solution [5, 26] and the negative diffraction, with front propagation [19]. It is 
natural to wonder which kind of solution exists in this limit. A linear stability study inform us that the two branch of 
the hysteresis cycle are linearly stable, contrary to the diffraction cases [20, 26]. Consequently, we expect to observe 
front solution, connecting the two branch of the bistable cycle. 

The numerical simulations from Eq. (1) have been performed for a zero cavity length, i.e. zero diffraction, reveal 
the presence of the two states, notwithstanding without front propagation [cf. Fig. 8(a)]. However the experiments 
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Figure 7: Temporal evolution of front position xo(t) corresponding to spatiotemporal diagrams of Fig. (4-6). (a) Experiment and 
(b) numerical simulations of the LL model, Eq. (1), used the same parameters considered in Eig. (4-6). Blue diamonds : location 
values extracted from the smoothed spatiotemporal diagrams. Dashed red curves: best fit using expression 14. Experimental 
fit parameters : a = 308 jum, b = 0.069, to = 35.2 s, numerical fit parameters : a = 306 /xm, b = 0.017, to = 107.7 s. 


realized for a zero diffraction cavity show two fronts, connecting the dark state with the bright state, with slow, 
asymmetric and irregular speeds [see Fig. 8(c)]. We conclude that the Lugiato-Lefever model is not enough to model 
the experimental observations for zero diffraction. Close to a = 0, the diffusion of the liquid crystal molecules cannot 
be neglected. Consequently, we complete the model by a second equation governing the spatial nonlinear refractive 
index evolution, the intracavity field E and the nonlinear refractive n are given by 


dE 

dt 


= Ein{x) — (1 + iA)E + inE — i \a\ 


d‘^E 

dr‘^ 


+ (x, t ), 


(17) 


d'^n 2 


(18) 


where a is the diffusion coefficient. The stochastic numerical simulations realized with this model show a good 
agreement with the experimental observations, as presented on Fig. 8. 


V. CONCLUSION 

We have investigated the formation of localized structures in the Lugiato-Lefever equation in a negative diffraction 
regime with an inhomogeneous injected beam. We first show how to generate an equivalent left-handed Kerr material 
in the visible range. Experimentally, we show that the nonlinear dynamical states appearing in a focusing Kerr 
Fabry-Perot cavity submitted to negative optical feedback are propagating fronts in inhomogeneous medium. 

We have performed a reduction of the Lugiato-Lefever equation to a simple bistable model with inhomogeneous 
injected beam. From this simple model we have derived a simple expression for the speed front. We have used the 
inhomogeneous spatial pumping in the form of Gaussian beam. The front moves and stops on a asymptotic position. 
The experimental trajectory of the front position under that forcing follows an hyperbolic tangent law that fully 
agrees with the prediction from a generic bistable imperfect pitchfork bifurcation model. 

Our analysis should be applicable to all fiber resonator [25] with an inhomogeneous injected beam. In this case the 
coupling is provided by dispersion. When dispersion and diffraction have a comparable influence, three dimensional 
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Figure 8: Spatiotemporal diagrams: (a) ID LL model, (b) ID LL model + liquid crystal diffusion, (c) Experiments. Experiments 
Iq = 455 W.cm~^, d = 0 mm, (f = —0.6 rad, Wx = 1400 /xm, Wy = 100 /xm, Ri = 81.8 %; R 2 = 81.4 %. Numerical simulation of 
LL model with Eq = 1.95, A = 3.0, a = 0, Wx = 1400 iim„ e = 0.4. 


localized structures can be generated [30, 31, 33-37]. These structures consist of regular 3D lattices or localized bright 
light bullet traveling at the group velocity of light in the material. We plan to extend our analysis to three dimensional 
cavities with an inhomogeneous injected beam. In this case, localized light bullet in the negative diffraction and in 
anomalous dispersion may be stable. 
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